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Abstract. We propose a geometric growth model for weighted scale-free networks, which is controlled by 
two tunable parameters. We derive exactly the main characteristics of the networks, which are partially 
determined by the parameters. Analytical results indicate that the resulting networks have power-law 
distributions of degree, strength, weight and betweenness, a scale-free behavior for degree correlations, 
logarithmic small average path length and diameter with network size. The obtained properties are in 
agreement with empirical data observed in many real-life networks, which shows that the presented model 
may provide valuable insight into the real systems. 
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1 Introduction 

Complex networks [ 1 |2I3| 1 [5] describe a number of real- 

■ life systems in nature and society, such as Internet [6], 
World Wide Web [7] , metabolic networks [8] , protein net- 

■ works in the cell [9], worldwide airport networks |10|llj . 
co-author networks [12113114115] and sexual networks [TB] • 
Since the publication of the pioneering papers by Watts 
and Strogatz on small-world networks [17j and Barabasi 
and Albert on scale- free networks [IS], modeling real- life 
systems has attracted an exceptional amount of attention 
within the physics community |1|2|3|4|5] . 

Up to now, the research on modeling real-life systems 
has been primarily focused on binary networks, i.e., edges 
among nodes are either present or absent, represented as 
binary states. The purely topological structure of binary 
networks, however, misses some important attributes of 
real- world networks. Actually, many real networked sys- 
tems exhibit a large heterogeneity in the capacity and the 
intensity of the connections, which is far beyond Boolean 
representation. Examples include strong and weak ties 
between individuals in social networks [12ll3ll4ll5j . the 
varying interactions of the predator-prey in food networks 
unequal traffic on the Internet [5] or of the passengers in 
airline networks |10|llj . These systems can be better de- 
scribed in terms of weighted networks, where the weight 
on the edge provides a natural way to take into account 
the connection strength. In the last few years, modeling 
real systems as weighted complex networks has attracted 
an exceptional amount of attention. 
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The first evolving weighted network model was pro- 
posed by Yook et al. (YJBT model) [20] , where the topol- 
ogy and weight are driven by only the network connection 
based on preferential attachment (PA) rule. In Ref. |21J, 
a generalized version of the YJBT model was presented, 
which incorporates a random scheme for weight assign- 
ments according to both the degree and the fitness of a 
node. In the YJBT model and its generalization, edge 
weights are randomly assigned when the edges are cre- 
ated, and remain fixed thereafter. These two models over- 
look the possible dynamical evolution of weights occur- 
ring when new nodes and edges enter the systems. On 
the other hand, the evolution and reinforcements of inter- 
actions is a common characteristic of real-life networks, 
such as airline networks [10111] and scientific collabora- 
tion networks [12|13|14|15j . To better mimic the reality, 
Barrat, Barthelemy, and Vespignani introduced a model 
(BBV) for the growth of weighted networks that couples 
the establishment of new edges and nodes and the weights' 
dynamical evolution 22 23J. The BBV model is based on 
a weight-driven dynamics [24] and on a weights' reinforce- 
ment mechanism, it is the first weighted network model 
19j;hat yields a scale-free behavior for the weight, strength, 
and degree distributions. Enlightened by BBV's remark- 
able work, various weighted network models have been 
proposed to simulate or explain the properties found in 
real systems |4l5l25l2ril27l28l29l3()] . 

While a lot of models for weighted networks have been 
presented, most of them are stochastic ;4|. Stochasticity 
present in previous models, while according with the ma- 
jor properties of real-life systems, makes it difficult to gain 
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a visual understanding of how do different nodes relate 
to each other forming complex weighted networks [31] . It 
would therefore of major theoretical interest to build de- 
terministic weighted network models. Deterministic net- 
work models allow one to compute analytically their fea- 
tures, which play a significant role, both in terms of ex- 
plicit results and as a guide to and a test of simulated and 
approximate methods |3ll32l33l34l35l36l37l38l3<)l40l41l42l43l44l45l46l47l4g] 
So far, the first and the only deterministic weighted net- 
work model has been proposed by Dorogovtsev and Mendes 
(DM) @H] . In the DM model, only the distributions of the 
edge weight, of node degree and of the node strength are 
computed, while other characteristics are omitted. 

In this paper, we introduce a deterministic model for 
weighted networks using a recursive construction. The model 
is controlled by two parameters. We present an exhaus- 
tive analysis of many properties of our model, and obtain 
the analytic solutions for most of the features, including 
degree distributions, strength distribution, weight distri- 
bution, betweenness distribution, degree correlations, av- 
erage path length, and diameter. The obtained statistical 
characteristics are equivalent with some random models 
(including BBV model). 




t=l 





t=2 



2 The model 



Fig. 1. Illustration of a deterministically growing network in 
the case of m = 2 and 5=1, showing the first three steps of 
growing process. The bare edges denote the edges of weight 1. 



The network, controlled by two parameters m and i5, is 
constructed in a recursive way. We denote the network 
after t steps by Q(t), t > (see Fig. [1]). Then the network 
at step t is constructed as follows. For t = 0, Q(0) is an 
edge with unit weight connecting two nodes. For t > 1, 
Q(t) is obtained from Q(t— 1). We add mw (m is positive 
integer) new nodes for either end of each edge with weight 
w, and connect each of mw new nodes to one end of the 
edge by new edges of unit weight; moreover, we increase 
weight w of the edge by mdw (5 is positive integer). In the 
special case 6 = 0, it becomes binary networks, where all 
edges are identical |35|53|42| . 

Let us consider the total number of nodes N t , the total 
number of edges E t and the total weight of all edges W t 
in Q(t). Denote n v (t) as the number of nodes created at 
step t. Note that the addition of each new node leads to 
only one new edge, so the number of edges generated at 
step t is n e (t) = n v (t). By construction, for t > 1, we have 

n v {t) = 2mW t -x, (1) 

E t =E t - 1 +n v (t), (2) 

and 

W t = W t -i(l + mS) +2mW t -i. (3) 

On the right-hand side of Eq. ([3]), the first item is the sum 
of weigh of the old edges, and the second term describe 
the total weigh of the new edges generated in step t. We 
can simplify Eq. §5§ to yield 

W t = (1 + mS + 2m)W t - 1 . (4) 

Considering the initial condition Wq = 1, we obtain 

W t = (1 + m5 + 2m)'. (5) 



Substituting Eq. ((5]) into Eq. {T]), the number of nodes 
created at step t (t > 1) is obtained to be 

n v (t) = 2m(l + m8 + 2m) t ~ 1 . (6) 

Then the total number of nodes present at step t is 
t 

= -^— 5 [(1 + m5 + 2m)* +6+1] . (7) 

Combining Eq. (J6]) with Eq. (f2|) and considering Eq = 1, 
it follows that 

Ek = ~~~ t [(1 + mS + 2mf + S] . (8) 
2 + o 

Thus for large t, the average node degree fc t = and 
average edge weight Wt — are approximately equal to 

2 and respectively. 

3 Network properties 

Below we will find that the tunable parameters m and 5 
control some relevant characteristics of the weighted net- 
work Q(t). We focus on the weight distribution, strength 
distribution, degree distribution, degree correlations, be- 
tweenness distribution, average path length, and diameter. 
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3.1 Weight distribution 

Note that all the edges emerging simultaneously have the 
same weight. Let w e (t) be the weight of edge e at step t. 
Then by construction, we can easily have 

w e (t) = (l + m5)w e (t-l). (9) 

If edge e enters the network at step r, then w e (r) = 1. 
Thus 

itfe(f) = (1 + mS)^ T . (10) 

Therefore, the weight spectrum of the network is discrete. 
It follows that the weight distribution is given by 

' fie(0) _ 5 + 2 



Here the first item accounts for the increase of weight of 
the old edges incident with i, which exist at step t — 1. The 
second term describe the total weigh of the new edges with 
unit weight that are generated at step t and connected to 
node i. 

From Eq. (| 1 5[) . we can derive following recursive rela- 
tion: 

Si(t) = (l + mS + m)si(t-l). (16) 



Using Si(tj) = 1, we obtain 

Si (t) = (1 + m<5 + m) t ~' i 



(17) 



E t 2[(1 + mS + 2mY + 5} 
P(w) = { n e (r) _ m(2 + 5) (1 + mS + 2m) T_1 



Et 







(l + mS+ 2m)' + 5 



and that the cumulative weight distribution |3l32j is 
-Pcum(w) 



for r = 0, 

for r > 1, 

otherwise 
(11) 



Since the strength of each node has been obtained explic- 
itly as in Eq. (fP7|) . we can get the strength distribution 
via its cumulative distribution [3|32j . i.e. 



E t 

fl<T 

(1 +mS + 2m) T + 5 



(1 + mS + 2m) 1 + 5 ' 
Substituting for r in this expression using t = t— '" 



(12) 



ln(l+m(5) 



gives 



ln(l + mi5 + 2m) 



(1 + m5 + 2myw i»(i+m») + 5 
(l + mS + 2m)' + <5 

ln(l + mi + 2m) 

u> in(i+ m «) for large t. (13) 



So the weight distribution follows a power law with the 
exponent lw = 1 + \ n(1+mi) ; . 



^cum(^) ^ ^ 



fX,<ti 

(1 + mS + 2m) u +5+1 
(1 + mS+ 2m)' + 6 + 1 ' 



(18) 



From Eq. ^TTJ>, we can derive U=t - ln(1+ '^ +m) . Substi- 
tuting the obtained result of ti into Eq. (I18p gives 

ln(l+mi+2m) 

(1 + m(5 + 2m)' s ^n+mi+m) +5 + 1 
(l + mS+ 2m)' +5 + 1 

lDtl + mi+im) 

s in<i+ m 6 +m) f or i arge ^, (19) 



-fcum( , 5) 



Thus, node strength distribution exhibits a power law be- 

ln(l+mg+2m) 
ln(l+m(5+m) 



havior with the exponent j s = 1 - 1 - in ( 1 + mS + 2m ) 



3.3 Degree distribution 



3.2 Strength distribution 

In a weighted network, a node strength is a natural genear- 
lization of its degree. The strength Sj of node i is defined 
as 

(14) 



Wl 3 ' 



where denotes the weight of the edge between nodes i 
and j, fli is the set of all the nearest neighbors of i. The 
strength distribution P{s) measures the probability that 
a randomly selected node has exactly strength s. 

Let Si{t) be the strength of node i at step t. If node i is 
added to the network at step ti, then Si(t;) = 1. Moreover, 
we introduce the quantity Asi(t), which is defined as the 
difference between Si(t) and Si(t— 1). By construction, we 
can easily obtain 



Asi(t) = Si (t) - Si (t - 1) 
= m5 Wij + m 



m5si (t — 1) + m Si (t — 1) . 



(15) 



The most important property of a node is the degree, 
which is defined as the number of edges incident with the 
node. Similar to strength, in our model, all simultaneously 
emerging nodes have the same degree. Let fcj(i) be the de- 
gree of node i at step t. If node i is added to the graph 
at step ti, then by construction ki(ti) — 1. After that, the 
degree fcj(i) evolves as 



ki(t) = ki(t - l)+msi(t - 1), 



(20) 



where msi(t — 1) is the degree increment Aki(t) of node i 
at step t. Substituting Eq. (fT7| into Eq. (|2"0|) , we have 

Aki(t) = m(l + m5 + mf- 1 -*' . (21) 

Then the degree fcj(t) of node i at time t is 



t 



h(t) = h(U) + > Aki(r]) 



ri=U+l 

(m5 + 1 + m)*--** + 5 



5+1 



(22) 
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Analogously to computation of cumulative strength dis- 
tribution, one can find the cumulative degree distribution 



ln(l + ma + 2m) 



(1 + mS + 2m)* [((5 + 1) k - S}~ ^ i+ m s+ m j 
(l + m6+ 2m)* +S+1 
6 + 1 



(1 + mS + 2m) t + 6+1 

for large t. (23) 

Thus, the degree distribution is scale-free with the same 

n(l+mi5+2ro) 
ln(l+m<5+m) 



ln(l + mi + 2m) 

[(6 + 1) k] li(l+«"5+m) 



exponent as 7 S , that is jk — 7s = 1 J — 



3.4 Betweenness distribution 

Betweenness of a node is the accumulated fraction of the 
total number of shortest paths going through the given 
node over all node pairs |50|13j . More precisely, the be- 
tweenness of a node i is 



bi = 



j=ti=tk 



a-jk 



(24) 



where (Xjk is the total number of shortest path between 
node j and fc, and <Jjk{i) is the number of shortest path 
running through node i. 

Since the considered network here is a tree, for each 
pair of nodes there is a unique shortest path between 
them [51152153] . Thus the betweenness of a node is sim- 
ply given by the number of distinct shortest paths passing 
through the node. From Eqs. (|2"T]) and (]2"2"]l . we can easily 
derive that for a < 9 the number of nodes with degree 
( mA '+i+™) + s which are direct children of a node with de- 

o + l 

gree +S is m(l + m6 + m) r ~ 1-Q . Then at time 

t, the betweenness of a ^-generation-old node v, which is 
created at step t — 8, denoted as bt{6) becomes 



b t {9)=C s t [N t -(C s t +l) 



e-i 



m(l + m6 + m) 1 



_ Y _ a (C^ + l 



(25) 



where Cf denotes the total number of descendants of node 
v at time t, where the descendants of a node are its chil- 
dren, its children's children, and so on. Note that the de- 
scendants of node v exclude v itself. The first term in 
Eq. ([25]) counts shortest paths from descendants of v to 
other vertices. The second term accounts for the shortest 
paths between descendants of v. The third term describes 
the shortest paths between descendants of v that do not 
pass through v. 

To find bt(6), it is necessary to explicitly determine the 
descendants Cf of node v, which is related to that of v's 
children via [53] 



E 



m(l + m6 + to)"- 1 {C T t - a + l) 



(26) 



Using C° = 0, we can solve Eq. ([26]) inductively, 
1 



r - 

Lt ~ 6 + 2 



[(mS + 1 + 2m) T - 1] . 



(27) 



Substituting the result of Eq. (g7]) and © into Eq. (g5 
we have 



bt(0) 



(m6 + 1 + 2m) 



t + T 



(£ + 2) 2 

Then the cumulative betweenness distribution is 
n v (//) 



(28) 



Pcum(b) 



_ (1+mS + 2m) t - T + 6 + 1 
_ (l + mS+ 2m) t + 6 + 1 
(l + m6 + 2m) t Nf 
(1 + mS + 2mY+T b ' 1 ' 

which shows that the betweenness distribution exhibits 
a power law behavior with exponent 7b = 2, the same 
scaling has been also obtained for the m = 1 case of the 
Barabasi-Albert (BA) model describing a random scale- 
free treelike network [51152] . 

3.5 Degree correlations 

Degree correlation is a particularly interesting subject in 
the field of network science [54 55 56 57 58 59], because it 
can give rise to some interesting network structure effects. 
An interesting quantity related to degree correlations is 
the average degree of the nearest neighbors for nodes with 
degree fc, denoted as k nn (k) [55156] , When k nn (k) increases 
with fc, it means that nodes have a tendency to connect 
to nodes with a similar or larger degree. In this case the 
network is defined as assortative [57 58 . In contrast, if 
fcnn(fc) is decreasing with fc, which implies that nodes of 
large degree are likely to have near neighbors with small 
degree, then the network is said to be disassortative. If 
correlations are absent, fc nn (fc) = const. 

We can exactly calculate fc nn for the networks using 
Eq. (]22[) to work out how many links are made at a par- 
ticular step to nodes with a particular degree. We place 
emphasis on the particular case of 6 — 0. Except for the 
initial two nodes generated at step 0, no nodes born in the 
same step, which have the same degree, will be linked to 
each other. All links to nodes with larger degree are made 
at the creation step, and then links to nodes with smaller 
degree are made at each subsequent steps. This results in 
the expression [34139] for fc = (m + l) t_ti 



fcnn(fc) 



1 



n v (ti)k(U,t) 

m-n v (t' t )k(t' t ,U - l)k(t[,t) 




t'i=t 



+ m-»H,(*i)fr(*i>*i » ( 30 ) 



=u+i 



Zhongzhi Zhang et al.: Recursive weighted treelike networks 



5 



where k(ti,t) represents the degree of a node at step i, 
which was generated at step t%. Here the first sum on the 
right-hand side accounts for the links made to nodes with 
larger degree (i.e. t\ < U) when the node was generated at 
ti . The second sum describes the links made to the current 
smallest degree nodes at each step t\ > ti. 

Substituting Eqs. © and ([12]) into Eq. ([30]) . after some 
algebraic manipulations, Eq. ([50)1 is simplified to 



2m + 1 
rn 

m + 1 



(m + 1) 



2lU 



2m + 1 
m 



m - 



1 



(t-U). 



(31) 



Thus after the initial step fc nn grows linearly with time. 

Writing Eq. ([3"T]) in terms of A:, it is straightforward to 
obtain 



2m + 1 



m 
m - 



(m+ l) 2 



1 



2m + 1 
m 



5m+l 



lnfc 



m m + 1 ln(m + 1) 



(32) 



t=0 



t=l 



(a) 



t=2 



t=3 



(b) 













, a. q? 




Q" 

















Fig. 2. (Color online) (a) The first four steps of binary network 
growth for the limiting case of m = 1 and 8 = are shown, 
(b) The network after t + 1 generation, Qt+i, can be obtained 
by joining three copies of generations t (i.e. Q[ , Qt , Qt ) at 
the two hub nodes of highest degree, denoted by A and B. 



Therefore, k nn (k) is approximately a power law function of 
k with negative exponent, which shows that the networks 
are disassortative. Note that k nn (k) of the Internet exhibit 
a similar power-law dependence on the degree k nn (k) ~ 
fc~", with to = 0.5 [55]. Additionally, one can easily check 
that for other values of 6 > 0, the networks will again be 
disassortative with respect to degree because of the lack 
of connections between nodes with the same degree. 



3.6 Average path length 

Most real-life systems are small- world, i.e., they have a 
logarithmic average path length (APL) with the number 
of their nodes. Here APL means the minimum number of 
edges connecting a pair of nodes, averaged over all node 
pairs. For general m and S, it is not easy to derive a closed 
formula for the average path length of Q(t). However, for 
the particular case of m = 1 and 5 = 0, the network has 
a self-similar structure, which allows one to calculate the 
APL analytically 

For simplicity, we denote the limiting network (m = 1 
and (5 = 0) after t generations by Qt- Then the average 
path length of Q t is defined to be: 



D, 



N t (N t - l)/2 



(33) 



In Eq. ([3"3"]) , D t denotes the sum of the total distances 
between two nodes over all pairs, that is 



D t = 



h3 ' 



(34) 



where dij is the shortest distance between node i and j. 

We can exactly calculate d t according to the self-similar 
network structure [44]. As shown in Fig. (2] the network 



Qt+i may be obtained by joining at the hubs (the most 
connected nodes) three copies of Qt, which we label Q t (Q) , 
a = 1,2,3 [42 60 . Then one can write the sum over all 
shortest paths D t+ i as 



Dt+! = 3A + At 



(35) 



where A t is the sum over all shortest paths whose end- 
points are not in the same Qt branch. The solution of 
Eq. ([31]) is 



D t =3 t - 1 D 1 +J23 t - r - 1 A T . 



(36) 



The paths that contribute to A t must all go through at 
least either of the two hubs (A and B) where the three 
different Q t branches are joined. Below we will derive the 
analytical expression for A t named the crossing paths, 
which is given by 



A t = A + At 



2,3 



tl,3 



(37) 



where Z\" ,/5 denotes the sum of all shortest paths with 
endpoints in Q^ and Qf '. 



If Q^ and Q[^ meet at an 
edge node, A"' p rules out the paths where either endpoint 
is that shared edge node. If Q^ and Q\ do not meet, 
A"' 13 excludes the paths where either endpoint is any edge 
node. 



By symmetry, A t ' 



Z\ 2 ' 3 , so that 



A t = 2A)' 2 + A\ 



.1,3 



where A]' 2 and A]' 3 are given by the sum 



A 1 ' 2 = 



E 



(38) 



(39) 
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and and 

A 1,3 = E ( 4 °) 

i^A, j+B 



'<■./ 



ieQr'.jeQf j#b 



A 1 ' 3 = E *• 

3 

E (^ + + 1 ) 



respectively. In order to find A\' 2 and A t ' 3 , we define 



4 ot 



ZeQl > di,A<di, Al , d 3:B >dj, Bl 



zeQf 3 zeQ?> i 6 Q«, jeQf, i^, j¥b 

<iz,A<<iz,B dz,A<dz,B di A>rfi,A 1 . dj B <dj,B 1 

d^= J2 dz - A > Ni t aT = E !> ( 41 ) + E (4i+%A + i) 

dz,A>dz,B d z ,A>dz,B d i A<di i A 1 , dj B <dj,B 1 

= 2(N t ~ l)(d? ear + 4 ai ) + ( N t ~ I) 2 , (45) 

where Z + 1 A. Since A and B are linked by one edge, for 

any node i in the network, di t A and <k,B can differ by at respectively. Having A^ 2 and Zi*> 3 in terms of the quan- 
most 1, then we can easily have d\ ot — <i" car + d[ al and tities in Eq. (|¥T|) . the next step is to explicitly determine 
N t = AT t near + N l t al + 1. By symmetry JV t near + 1 = iV 4 far . these quantities unresolved. 

Thus, by construction, we obtain Considering the self-similar structure of the graph, we 

can easily know that at time t + 1, the quantities d' t + ar and 



N t = 2(iV t noar + 1). (42) 



are related to each other, both of which evolve as 

7near jfar , r> jncar 



J ( f t l ; a 1 r = 4 ar + 2rf t loar , 

Combining this with Eq. ([7]), we obtain partial quantities i ^far _ ^ncar , ^yfar \^) 

in Eq. gj) as * * ' 



From the two recursive equations we can obtain 
- 1 = AT" = J (3* - 1) ■ (43) f , ncar 1 



Now we return to the quantity A\ 2 and A]' 3 , both of ^ 12 

which can be further decomposed into the sum of four 

terms as Substituting the obtained expressions in Eqs. (|43|) and 



^near = _ £_g + gl+t + ^ . gt J 
I 12 

^ ar = _ (3 + 3 2+t + 2t • 3*) . 



(47) 



(|47j) into Eqs. l(34)). (|45j) and [|38]). the crossing paths A 
12 ^ s obtained to be 

A' = 2^ d *j 

ieQf'.ieQf A = 7-9* + 2t-9*. (48) 



E (^m + 



Inserting Eq. (|48j) into Eq. ((36)) and using D\ — 10, we 
have 

fQ?K 3&Q?\ idfA D = 3-1+t (i + 2 . 3* + 1 • 3*) . (49) 

<ii,A><ii,A 1 , dj,A>dj,B V ' V ' 

Substituting Eqs. and (jl9")l into (p?3")) , the exact expres- 
sion for the average path length is obtained to be 



E ( di < A + d 3A) 

A<di, Al , di.AXij.B _ 2(1 + 2-3* + t- 3* 



E Ka+M "* 3(1 + 3') ■ (50) 



In the infinite network size limit (t — » oo), 



di,A><ii,A 1 , dj,A<dj t B 

+ „ S ( ^ + ^' + (51) 

di } A<di } A 1 , dj,A<dj t B 

ncar , which means that the average path length shows a loga- 

— 2(JVt — l)(d t + d t ) , (44) r ithmic scaling with the size of the network. 
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3.7 Diameter 

Although we do not give a closed formula of APL of Q(t) 
for general m and 5 in the previous subsection, here we will 
provide the exact result of the diameter of Q(t) denoted 
by Diam(Q(t)) for all m and 8, which is defined as the 
maximum of the shortest distances between all pairs of 
nodes. Small diameter is consistent with the concept of 
small-world. The obtained diameter scales logarithmically 
with the network size. Now we present the computation 
details as follows. 

Clearly, at step t = 0, Diam(Q(0)) is equal to 1. 
At each step t > 1, we call newly-created nodes at this 
step active nodes. Since all active nodes are attached to 
those nodes existing in Q(t — 1), so one can easily see that 
the maximum distance between arbitrary active node and 
those nodes in Q(t— 1) is not more than Diam(Q(t— 1))+1 
and that the maximum distance between any pair of ac- 
tive nodes is at most Diam(Q(t — 1)) + 2. Thus, at any 
step, the diameter of the network increases by 2 at most. 
Then we get 2(t + 1) as the diameter of Q(t). Note that 
the logarithm of the size of Q(t) is approximately equal to 
tln(l + mS + 2m) in the limit of large t. Thus the diame- 
ter is small, which grows logarithmically with the network 
size. 



4 Conclusion 

In summary, we have introduced and investigated a deter- 
ministic weighted network model in a recursive fashion, 
which couples dynamical evolution of weight with topo- 
logical network growth. In the process of network growth, 
edges with large weight gain more new links, which oc- 
curs in many real-life networks, such as scientific collab- 
oration networks |49ll2ll3ll4ll5j . We have obtained the 
exact results for the major properties of our model, and 
shown that it can reproduce many features found in real 
weighted networks as the famous BBV model [22123] . Our 
model can provide a visual and intuitional scenario for the 
shaping of weighted networks. We believe that our study 
could be useful in the understanding and modeling of real- 
world networks. 
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